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1. Introduction and definitions 

The Harish Chandra-Itzykson-Zuber integral P],EJ plays a key role in matrix models 
0,|3],f|.|5|,[| , in particular in relation to (discretized) two-dimensional quantum gravity. It 
also displays interesting connections with free probability theory PJToH. It is defined 
by 

I N (A,B)= d0e iVtrAOBOt 



JtlEU(N) 

where A and B are N x N matrices and dO is the normalized Haar measure on U(N). 
There are many ways to evaluate this integral |I|P,^, pT| . |I^ , |13| ; the result depends only on 
the eigenvalues Gtj of A and 6j of B, and is given by 

A(a i )A(6 i ) 

Here A(-) denotes the Vandermonde determinant: A(a^) = Yli<j( a i~ a j)- ^ n au equations, 
constants that only depend on A" will be ignored; they can be restored when required. 

For practical applications it is important to look at the large A" limit of the HCIZ 
integral: it corresponds to the "planar" limit of matrix models. One considers a sequence of 
N x N matrices A and B, whose spectral distributions converge to some fixed distributions 
as A" goes to infinity. It is known (and in fact, proved rigorously ||14|| ) that the following 
quantity, the "large A" free energy" , is well-defined: 

and depends on the spectral distributions of A and B only through their moments: 

e„ = lim —A q L = lim —B q q > 1 (1.4) 

We also define for future use the generating series of these moments, that is the resolvents 
of A and B: 

G(a)= lim — — — = l/a + 6Ja q+1 (1.5a) 
G(b) = lim | -J- = l/b + J £e q /W+ 1 (1.56) 

As a function of the 9 q and the 6 q , F has a power series expansion around 0. We shall 
investigate here the determination of the coefficients of this expansion, as well as the more 
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general question of the evolution of F as one varies the 6 q and 9 q . We shall see that this 
evolution is governed by an infinite set of partial differential equations, a particular scaling 
of the Toda lattice hierarchy. 

More precisely, we shall rewrite in section 2 the HCIZ integral as a unitary two- 
matrix model, then show in section 3 that the latter satisfies the Toda lattice hierarchy; 
in section 4, we shall take the large N limit which relates F to the dispersionless Toda 
lattice hierarchy; finally, we shall apply this technology to the expansion of F in section 5. 
Section 6 is devoted to a summary of results. 



2. HCIZ integral as a unitary two- matrix model 

We shall first use a trick to rewrite the HCIZ integral as an integral over two matrices. 
The derivation can be performed either by computing residues in the two-matrix integral, 
or by character expansion; we choose the latter. We thus start with the following identity, 
valid for any N x N matrix M: 



e iVtrM 



J2cRd RX R(M) (2.1) 



R>0 



Let us explain the symbols used. R stands for an analytic irrep of GL(N) with positive 
highest weights, or the associated Young diagram; more explicitly, it can be described by 
its "shifted" highest weights hi, 1 < i < N which form a strictly decreasing sequence (they 
are related to the usual highest weights mi by hi = m 8 + N — i). xk(M) is the character 
associated to the irrep R and the matrix M. cr and d R are coefficient of the expansion: 

cr = ^ N d R = A(hi) (2.2) 

where \R\ = £V hi — N(N — l)/2 is the number of boxes of R. These explicit expressions 
of the coefficients play no role in what follows since the trick we use is quite general.0 

Applying the identity ( |2.1| ) to the definition of the HCIZ integral (|1 . 1| ) and using 
orthogonality relations for matrix elements of irreps yields 

I N (A, B) = J2 c rX r(A)xr(B) (2.3) 

R>0 



1 One can treat similarly any expression of the form a(V-)A(I> J -) > c ^ -^l' (H3)- However here we 



do not use Eq (L2) - in fact, we prove it along the way. 
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As a side remark, we note that using Weyl's formula for characters (xii(A) = det(a^ J ) / A(a^)) . 
recombining the two resulting determinants into one single determinant and performing 
the summation, one recovers Eq. ( |1 . 2| ) . 

We now want to compare expression ( |2.3| ) with the partition function of a matrix 
model of two N x N matrices U and V: 

r N \t„ i q ] = ffdrnVe^*^ + E„>i «>V + MJ-»V-' p 4) 

Here the t q are as yet unknown coefficients; they must however satisfy sup g l^l 1 / 9 < m < oo 
and sup 9 l^l 1 ^ 9 < m < oo; this renders the summation convergent provided that the 
spectral radii of U and V are less or equal to 1/m, 1/m. Since the integrand is analytic in 
U and V the exact contours of integration are irrelevant; for the sake of definiteness, we 
choose U and V to be unitary up to a multiplicative constant: UU^ = m -2 , VV^ = fh~ 2 . 
The measure of integration is simply the normalized Haar measures for m U and m V. All 
these details can be safely ignored if one remembers that in the integral (|2.4| ) only the 
contribution coming from poles at the origin should be picked up. 

Note that we have defined ( |2.4[ ) so that N does not appear anywhere explicitly in the 
action. It is however convenient to rescale the matrices to restore the usual N dependence 
of matrix models: 

t*K i,\ = //dUdVeS„>. t,N-'HrW + E„>i t,N-IH^ + WtrU-'V' (2 5) 
Next we apply one of the forms of the Cauchy identity: 

^ ,/2trU * = J2 SR(t,N-«<*) XR (lJ) (2.6) 

R>0 

where sr is the Schur function associated to the Young diagram R. 

Plugging twice (|2.6|) as well as the identity (|2.1|) into (|2.5| ) results in: 



tnKQ = s Rl {t q N-il 2 )s R2 {i q N-«/ 2 )c R d R U dUdVxRA^XR^XRiU^V- 1 ) 

R,R 1 ,R 2 >0 J J 

(2.7) 

Again, orthogonality relations imply that 

r N [t q ,t q ]= J2 CRS R {t q N-^ 2 )s R {i q N-^ 2 ) (2.8) 

R>0 
#rows<JV 
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Note that we have written explicitly the constraint that the Young diagram R should have 
no more than iV rows, since contrary to prior expressions the summand in Eq. (|2.8| ) is not 
necessarily zero if the number of rows exceeds N. 

Comparing Eqs. ( |2.8[ ) and ( J2 . 3| j , we see that the two expressions are equal (up to a 
numerical factor) iff 

t q N~ q / 2 = -trA 9 t q N- q / 2 = -trB 9 (2.9) 
q q 

or in terms of the normalized moments introduced before 

t =N i/2+ik l (2.10) 
q q 

The constraint on the number of rows of R is automatically satisfied if the t q are of the 
form (|2.9| ). However, the expression ( |2.8| ) has the advantage that it is defined for arbitrary 
(independent) t q at finite N, which is not the case of the original HCIZ integral (where 
there are only iV independent t q ). Also note that the condition of finiteness of sup 9>1 l^q] 1 / 9 
and sup g>1 l^l 1 / 17 is automatically satisfied at finite N if t q is of the form ( |2.9|) , and remains 
true as N — > oo provided the spectra of A and B remain bounded. 



3. Biorthogonal polynomials and Toda lattice hierarchy 

We have found in the previous section that calculating the HCIZ integral is equivalent 
to solving the matrix model given by Eq. (|2.4|). A convenient way to do so, at least 



formally, is via biorthogonal polynomials. This is standard material, and is very similar to 
derivations found in the literature ||15|| . tn will then be identified with a tau function of 
the 2D Toda lattice hierarchy. 

3.1. Setup 

We start with the following determinant form of the partition function (|2.4| ) , obtained 
by diagonalizing U and V: 

TN = det {??2^2™ UVe q - q ~ ) (3 - 1} 

\J J zntuzntv / Q<i,j<N-l 

where once again the contours of integrations are sufficiently small circles around the origin. 
This suggests to introduce a non-degenerate bilinear form on the space of polynomials 

by 

la \p)=<P<P T^—T^—pMqMe^^ 1 q q (3.2) 

X 1 ' J J 27TIU ZTTIV 



so that tn = detM N , with M N = (m^), < i,j < N — 1 and m^- = (i^J-u- 7 ). 

Next we define normalized biorthogonal polynomials q n (v) and p n (u) with respect to 
the bilinear form above on C[v] x C[u], that is (q n \Pm) = <^nm- We write: 

n n 
Pn(«) = ^PknU k q n (v) = ^2qknV k (3.3) 



k=0 k=0 



-1 
n • 



There is an arbitrariness in the normalization which is fixed by assuming p nn = q nn = h 

Define upper triangular matrices Pn = (Pkn) and Qn = {qkn)i < n, k < N — 1; and 
Pn = Qn 1 ■ Orthonormality of these polynomials is equivalent to writing 

Q T N M N P N = 1 (3.4) 

or Mn = PnPn , with P^ 1 upper triangular and Pn lower triangular. In particular 
r N = detM N = UtL~o 1 hl 

Next, introduce the semi-infinite matrices M = M^, P = P-^, Q = Qoo, P = -Poo 
and the shift matrix Z: Z^ = j+i • Set 

U = P~ 1 ZP V = Q~ 1 ZQ (3.5) 

U (resp. V) is the matrix of multiplication by u (resp. v) in the basis (p n (u)) of C[u] (resp. 
(g n («)) of C[v]). Also set £/ = V T = p- 1 Z T P. By definition of the rriij, we have 

f=M Z - ^=Z«M (3.6) 

We then easily find, using M = PP -1 , 

_p-i«P + p-i^ = U » _p- 1 E + p-'E = f /» (3.7) 

P~ 1 W~ * s u PP er triangular, whereas P -1 ^ is lower triangular, and they have opposite 
diagonal elements (and similarly for derivatives with respect to t q ); so that if one defines 
symbols (• • •)+ = (• • -) >0 + \{- ■ -)o and (• • •)_ = (• ■ -) <0 + §(• • -)o for lower and upper 
diagonal parts respectively (it is convenient to include one half of the diagonal part, though 
this does not make them projections), we have 

dP dP 

wr -P(u<). Wq =-r(u')- (3*0 
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and finally 



dU dll 

— = -[(U9) + ,U] ^- = [(U*)-,U] (3.9a) 

dU ~ dU 

0f q =-[(U 9 )+,U] ^ = [([/«)_,£/] (3.96) 



Equations (|3.9|) are equivalent to equations of Zakharov-Shabat type [16 



^ ur )+ - w r {uq)+ + (tr)+] = (3 ' 10a) 

J-(^) + + j^{V q )- - i(U q )-, (U r ) + ] = (3.106) 
J-(tn_ - §- r (U q )- ~ i(U q )-, (U1-] = (3.10c) 



which ensure compatibility of Eqs. (|3.8|) . If we expand U and U in the following way: 

oo 

U = rZ + J2 u kZ Tk (3.11a) 

oo 

U = Z T r + Y,Z k v k (3.116) 

k=0 

where r = diag(r n ), the u k = diag(«fc ;n ) and the v k = diag(f/ c; n) are diagonal matrices, 
then Eqs. flOoD are the (semi-infinite) Toda lattice hierarchy of differential equations for 
the coefficients of U and U. 

3.2. Examples 

The simplest equation, the usual Toda lattice equation, which is the case q = r = 1 in 
Eq. ( |3.106|) , can be obtained directly from the determinant form (|3.1] ) by applying Jacobi's 
determinant identity, which yields T n+1 r n -i = r n d ^ Q ^ r n - ^r n -^r n , or 

= mM loeT " (3 ' 12) 

The left hand side is nothing but (h n /h n -i) 2 = r 2 . This implies 

r 2 n+1 +r 2 n _ 1 -2r 2 n = -^logr^ (3.13) 

Ot\Oti 
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which is the Toda lattice equation in terms of the log r\ . In the next section we shall see 
how to use this equation in the large N limit. 

Another equation one can derive (case q = 2, r = 3 of Eq. (|3.10a| )) is the usual 

d 2 
i>t- 



o2 

Kadomtsev-Petviashvili (KP) equation, which is given in terms of x — 2^W logrjv: 



3.3. String equation 



So far we have not used the specific form of the interaction between our two matrices 
in ( |2.4[ ). This occurs when one tries to determine the string equation, that is an additional 
relation which fixes a particular solution of the Toda hierarchy. Here it is very easy to 
derive the string equation by arguments similar to those used in the standard two-matrix 
model. We write that 



, , du dv d 

= <t> <t> u — 

2niu 2niv du 



Pj(u)qi(v)e^ 



(3.15) 



expand and obtain 

UD + J2 at q uq -U^U' 1 = (3.16) 
«>i 

where D is the derivative with respect to u in the basis p n (u), and U~ x is a particular 
left inverse of U whose matrix elements are defined by Eq. (|3.15|) , and similarly U~ x is a 
particular right inverse of U. Using [D, U] = 1 and U~ 1 U = 1, we derive from Eq. ( 3. 16] ) 
the string equation: 

[U-\U- 1 ] = 1 (3.17) 

Let us note that Eq. ( |3.17|) is not the usual string equation of the hermitean two-matrix 
model. 



4. Large iV limit 

4-1. Classical/ dispersionless limit 

Next we want to consider the large iV limit. As always in matrix models, this identifies 
with a classical limit for the operators involved, the small parameter being h = 1/N. 
Therefore in this limit, the various operators U, U, Z become classical variables u, v, z, 
up to some rescalings which will be detailed later. 
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What happens to the Toda equation? The large N limit corresponds to the disper- 
sionless limit of the Toda hierarchy [T7],|18]]. In this limit, operators of discrete difference 



in the index n (e.g. Eq. ( |3.13| )) can be replaced with a differential operator djdv with 
v = n/N. Furthermore, the logarithm of the tau-function tn is required to be of order iV 2 
(cf our Eq. Q1.3Q ), with a smooth 1/N expansion,! which leads to some simplifications in 
the differential equations. We shall only write down a few examples of these dispersionless 
equations; we refer the reader to e.g. fT7j , |18l , |2~0|1 for a complete and more rigorous descrip- 
tion of the dispersionless hierarchy. Let us simply mention that, as mentioned above, we 
are dealing with a classical limit, that is the commutators become of order h = 1/N, and 
their leading term defines a Poisson bracket. Equations ( |3.9| )-( j3T7|) can be rewritten in 



the large iV limit, once rescaled, by simply replacing commutators with Poisson brackets. 
In particular, the string equation ( |3.17| ) becomes: 

{u-\v- 1 } = l (4.1) 

(where now w _1 and v _1 are ordinary inverses since in the large iV limit all operators are 
invertible) . It is worth noting that the same string equation occurs in the interior Dirichlet 
boundary problem [|21[| . Also, it is identical to the string equation of in which the 
radius has the unphysical value R = — 1. 

4-2. Scaling. Some examples 

In order to give some examples, we need to give some additional details on the par- 
ticular scaling of the variables that is appropriate for our model. First, we have for our 
operators U, U, Z: (cf the rescaling from Eq. (|2.4|) to Eq. (|2.5|) ): 

U ~ N~ 1/2 u U ~ N~ 1/2 v Z~z (4.2) 

where u, v , z have a finite limit when iV — > oo and are related by: 

oo 

u = rz + J ^2u k z~ k (4.3a) 

oo 

■ 1 + ^2v k z k (4.36) 



v = r z 

k=0 



The assumption of smoothness is less innocent than it seems; it is the equivalent of the "single 



cut" hypothesis in matrix models (cf a similar remark in [19| the context of ID Toda). However 



in formal power series around t q = t q = this hypothesis is certainly valid. 
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in terms of appropriately rescaled coefficients (in particular ~ riV -1 / 2 ). We have 
already derived the rescaling of the moments, cf Eq. ( |2.f 0| ). 

In what follows, we shall freely change from original to rescaled variables, using the 
same letters by abuse of notation. In most circumstances this does not cause any confusion. 

For example, let us see what happens to the KP equation ( |3.14| ). In the large N limit, 
performing all the rescalings one term drops out and we find 

2S + ^f-2#+x#l=0 (4.4) 



del V de 3 00! 

which is just the dispersionless KP equation (or Khokhlov-Zabolotskaia equation). 

Let now see what happens to the Toda lattice equation (|3.12| ) in the large iV limit. 



Here, it is important to notice that an extra ingredient is required if one considers equations 
in which variations with respect to n appear. It is our scaling hypothesis, based on the 
fact that we consider only the "planar" large iV limit (first term in the 1/N expansion) 
and on the scaling of the t q and i q . Considering Eq. (|1.3| ), we write that: 

Cn e n 2 F[n-^ 2 -V ? , n-i/^qt,} + 0(1) (45) 



The constant c n can be determined by setting t q = t q = in Eq. (|3Tl"|); one easily finds 

More explicitly, if one introduces the dispersionless tau function F[i>,9 q ,9 q ] = 
limAr^ 00 (-^- logr n + \v 2 logiV — |) with v = n/N, then we have the scaling 

F[u,e q ,e q ) = -\u 2 \o g u + - A {u 2 - 1) + ^F[V" 9/a_1 , V" 9/a_1 ] ( 4 -6) 

where our old function F[9 q , 9 q ] is recovered by setting v = 1. 
We now rewrite the Toda equation ( |3.12|) : 

e d^* = (4.7) 

which is just the dispersionless Toda equation (integrated twice); then apply the scaling 
(fO|). The final result is the following partial differential equation for F: 

exp (2F + + l)(f - 2)(e q d q F + 9 q d q F) 

+ (f + !)(§ + l)(e q e r d q d r F + 29 q 9 r d q d r F + 9 q 9 r d q d r F)^j = dxdxF (4.8) 

<j,r>l 
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with d q = d/d9 q and d q = d/d9 q . This is in general a fairly complicated equation. 
Application. Assume that all 9 q and 9 q are zero except 9\ and 9\. This greatly simplifies 
Eq. ( |4.8|) . Note furthermore that by obvious homogeneity property, F only depends on the 
product x = 9\9\. Rewriting Eq. ( |4.8|) for F(x) leads to an ordinary differential equation: 



e 2F + 9x*F" =F , + xF „ (49) 

with the initial condition -F(O) = 0, which can be solved: if £ = xdF/dx, then £ is the 
solution of the third degree equation 

16£ 3 + 8£ 2 + (1 - 36x)£ + x(27x - 1) = (4.10) 

that vanishes at x = 0. We have the following expansion: 

?= fv +1 (y (4 . n) 

s ^ 2n+l)! v ' 

n=0 

Curiously, the number ( n+ ^7(2n+i)! a ^ so a PP ears m the context of the enumeration of 
planar maps;! it counts trivalent rooted maps with 2n nodes ||23|| . This suggests that a 
direct combinatorial proof of this result might be possible. As a corollary, the HCIZ integral 
with #i, 9i ^ is in the universality class of pure 2D quantum gravity: the singularity of 
£(x) closest to the origin is of the form £ = (x — x c ) 3 / 2 + ■ • • (with x c = 2/27, £ c = 1/12). 

More generally, if one turns on a finite number of 9 q and 9 q , one can show using the 
formalism developed in this section that F can be expressed in terms of the solution of an 
algebraic equation. This is because the string equation implies that a = 1/u and b = 1/v 
are polynomials in z , z: 

I+i e+i 

q=l q=l 

where £ = max-f^g 7^ 0}, I = ma,x{q\9 q 7^ 0}, and as will be shown explicitly in the 
sections below (cf Eq. (4.18)), the knowledge of a and b allows to compute derivatives of 
the free energy .@ 



3 The author would like to thank J.-B. Zuber for pointing this out to him. 

4 Note that since a, b have a rational parameterization, the curve (a, b) has genus zero. 
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For example, if only 6>i, #2, • • •, 9g, 6\ are non-zero, by homogeneity one can set 6>i = 1; 
then, x = 24a? = 2?/> 4 , and r 2 = -f ^ = °0 are given by the degree 2£ + 1 equation 

-1 + ^ + Y {-l) q ^r6 q ip 2q+1 = -1 + V - 2#iV> 3 + 60 2 V> 5 - 20# 3 ^ 7 + • • • = (4.13) 

q~l ( ql > 

with the initial condition -0 = 1 for all # 9 zero. One can check that if I > 3, x = 2-0 4 does 
satisfy the partial differential equation (fOj). 

It is expected that by turning on an arbitrary number of 6 q and 9 q , one should be 
able to recover the critical behavior (i.e. singularities of F and its derivatives) of all c < 1 
minimal models coupled to gravity. 

4-3. Derivatives of F with respect to the t q 

Our goal is now to express the first few derivatives of F with respect to the t q and 
the t q . Let us define the differential operator 

q>l q 

In the bosonic language, this is the positive modes of the full bosonic field (f>(a) = 
J2 q >i + loga — J2 q >i tq a ~ q ■ One could define similarly V(6). 

Our first task is to compute V(a)F (see also [p4|,|7|| for a "saddle point" ap- 
proach to this question). We note that in terms of the matrix model ( |2.4j ), V(a)F = 
- limAr^oo (trlog(l - aU)), or 

-^-V(a)F = u 2 D(u) - u (4.15) 

where D(u) = limjv^oo ^tr - ^ y is the resolvent of U, and with the identification u = 1/a 
(we recall that appropriate scaling as N — > oo is implied). Next we use the following 
standard identity: 

p N (u) = h- 1 (det(u~V)) (4.16) 



where the average is again over N x N matrices with the measure given by ( |2.4j) . This 
expression of pn(u) can be proved by checking that (q n \pN) = o~nN using it. Considering 
the logarithmic derivative with respect to u of Eq. (|4.16[ ) in the large iV limit, we see that 



D(u) is simply the operation of derivative with respect to u = 1/a in the basis of the p n (u). 
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This operator was already introduced in the discussion of the string equation; Eq. ( |3.16| ) 
becomes, in the large iV limit, 

D(u) =u- 2 (v~ 1 -G(u) +u) (4.17) 

where we recall that G(u) = u + Yl q >i Qt q u q+1 is the resolvent of the original matrix A. 
Combining Eqs. ( |4.15| ) and ( |4. 17| ) leads to 

4-V(a)F = b(a) - G(a) (4.18) 
da 



where b= 1/v, a = 1/u and u and v related by Eqs. Q13| ), i.e. b(a) is obtained by composing 
the series b(z) = l/v{z) and the inverse series z(a) of a(z). Note that b(a) —G(a) is simply 
the non-negative part of the Laurent expansion of b(a) = 4-(f)(a) in powers of a. 

Note that to derive the identity ( fll8| ), we had to use the string equation under 



the form ( |3.16|) . In contrast, the second derivative V(ai)V(a2)-F is "universal" i.e. in 
our context it means that it is true for any solution of the dispersionless To da hierarchy 
regardless of the string equation. Derivation of V(ai)V(a2)-F requires two identities. We 
shall prove the first one, which is an expression for z. We note that since Z is the shift 
operator, one should have in the large N limit z = P n +i(u)/P n (u), that is taking the 
logarithm and using Eq. (|4.16| ), 



logz = logM-logr- V — F (4.19) 

q ovdt q 

q>l 



where it is recalled that the dependence on v of F is given by Eq. 

The second relation we use the dispersionless limit of the Fay identity [p!8[ , which we 
shall not prove here and write under the form (see ||25|| ): 



y, U ± U 2 Q* F y, -q p y, -q Q* „ 

f Ul _ U2 \ e ^Q,r>l q r dt q dt r r _ UlQ 2^q>l a l dvdt q r _ ^ g 2^q>l a 2 dvdt q r 

(4.20) 

Combining Eqs. QL19|) and ( [1.20|) , using the variables a\ = 1/wi, a2 = l/«2 and the power 
series inverse z(a) of a(z) leads to 

V(a0 V(a 2 )F = log ~ ^ + logr (4.21) 

1/ai - l/a 2 

(cf also a similar formula in 0). 

One can prove a formula for mixed derivatives: V(a)V(6)F = — log(l — z(b)/z(a)), 
but it will not be needed here. 

Formulae for higher numbers of derivatives become more and more complicated, (see 
in particular "residue formulae" for three derivatives |2C[j), and we shall not write them 
down here. 
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5. Expansion around of the HCIZ integral 

We now come to the main goal of this paper: to apply the equations above to the 
expansion of the large N HCIZ free energy F[6 q , 6 q ] in powers of 9 q and 9 q . We shall 
consider the problem in an asymmetric way: we shall apply the operator V(a) defined in 
the previous section to the free energy, that is consider derivatives with respect to the 9 q , 
and then set 9 q = and see how the result depends on the 9 q . The identities obtained 
above involve the functions b(a) = l/v(a = 1/u) and z(a); let us see how to compute these 
functions when all t q are zero. 

If we set t q = in def. (|3.2| ) of the bilinear form, we see that no positive powers of u 
are generated, so that rriij = (u 4 it- 7 ) satisfies 

m ^ = H^ 3>i (5 - 1} 



where we have defined s q by (p(v) = exp ( J2 q>1 t q v q ) = J2 q > s q v q ; they are Schur func- 
tions corresponding to single row Young diagrams. 

Let us first use the fact that M is upper triangular. Since M = PP~ l with P upper 
triangular and P lower triangular, P is diagonal i.e. q n (v) = Vn\v n . The normalization 
h n = (n!) -1 / 2 implies that r n = n -1 / 2 , that is after rescaling r = 1. Finally, we conclude 
from the definition (|3.5|) of V that all the coefficients v k in its expansion (or of its transpose 
U, Eq. (3.11)) are zero, or in the large iV limit 

v = z~ x (5.2) 

(note that this is consistent with Eq. ( |4.12| )). We have now reduced the problem to a single 
function b(a) = l/v(a) = z(a); its calculation was performed in using a saddle point 
method, but it is particularly simple to rederive it in our framework. We start again from 
Eq. (|5.1| ) and notice that one can easily invert this matrix, which yields the coefficients of 
P = M~ 1 P '. If we introduce the e q defined by ip(y ) _1 = exp ( — J2 q >i tq vq ) = J2 q >o &q v<1 
(they are (— l) q times the Schur functions corresponding to one-column Young diagrams), 
then expansion of the identity ( / 9(f)(/?(f) _1 = 1 leads to 

f j < i 
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We now define the matrix A = (a^) by A = P 1 Z T P 1 which according to defs. (|3] 
is a left inverse of U. We have 




Vi-<0' + 1 ) !g fc-i- 



kl * — ' j 

j,i<j<k-i J 



■k-l-i 

= \lbiYl + « + l)s 9 e fc _i_ 
9 =o 



(5.4) 



This time we recognize the expansion of ip'(v)ip(v) 1 = J2 n >o vH 1 E^o^^"-?' but we 
also have ip' (v^iy) -1 = J2 n >i nt n v n ~ 1 , so that 



k < i + 1 

\fk _ k = i + 1 

(k — i — fe > i + 1 



(5.5) 



This is an exact expression at finite N. In the large iV limit this can be simply rewritten, 
performing the appropriate rescalings and remembering that z = b, 



a 



9>1 



(5.6) 



We recognize the resolvent of the original matrix B (Eq. (1.5)): a(b) = G(b) (again this is 
consistent with the "dual" equation of ( pTTgp , ^V(6)F = a(b)-G(b)). Finally b{ a = 1/u) is 
the functional inverse of a(b); it is well-known that 6(a) is nothing but the generating series 
of free cumulants fh q of B [|10| (cf also the expression of connected correlation functions in 
matrix models, as in Eq. (31) of 



(5.7a) 
(5.76) 



j oo 

6(a) = - + y^m q a q ~ 1 



q=l 



v (q + J2 i ai-2)l T[ (-d 

^ (q-iy. J-- 1 - 



We can now apply the formulae obtained in previous section. We recall that V(a) = 



V >i — First we rewrite Eq. (|EI8D for = 



oo 



— V(a)F = 6(a) = > ///,// 

da a ^— ' 

9=1 



9-1 



(5.* 
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This equation can be found in 0, but it can already be extracted from 0; see also [|j for 
a combinatorial proof. 

The second derivatives are given by Eq. ( |4.21| ) with z{a) = 6(a): 

6(ai) - 6(a 2 ) 



V(ai)V(a 2 )F = log 



1/ai - l/a 2 



9-1 



lo § ( 1 - m i a i a 

<J=2 k=l 



k^q — k 
2 



(5.9) 



Since the expression is in fact known for all t q , one can differentiate once more and obtain 
the expression for the third derivatives at t q = with little extra work. First one computes 
the expression of b(z) at first order in t q via Eq. (|3.96|) ; one obtains after inversion 



r z 



q>l k>0 







(5.10) 



where the coefficients c q k are best expressed using their generating function 

> 9 d6(a) 



— C qk 



9>1 



da 



-6(a) 



-fc-2 



(5.11) 



Secondly one uses Eqs. ( |4.18| ) and (|4.21[ ) to deduce 



V(a 2 )6( ai ; 



d 



log(«(oi) - z{a 2 )) 



(5.12) 



Acting with V(a3) on Eq. ( [4.21|) , setting t q = 0, using identities ( |5.10|) - (|5.12|) and per- 
forming simple manipulations results in 



V(a 1 )V(a 2 )V(a 3 )F 



d6(ai)/dai 



+ 



d6(a 2 )/da 2 



+ 



(6(ai) - 6(a 2 ))(6( ai ) - 6(a 3 )) (6(a 2 ) - 6( ai ))(6(a 2 ) - 6(a 3 )) 
d6(a 3 )/da 3 



(6(a 3 )-6(a 1 ))(6(a 3 )-6(a 2 )) 



+ 1 



(5.13) 



One can go to higher derivatives by further differentiation, but most likely, the expressions 
become more and more cumbersome. 



6. Summary of results 

In this paper we have investigated the expansion of the large iV HCIZ integral free 
energy in terms of powers of the moments of the matrices. We have shown that this free 
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energy is the dispersionless tau function for the dispersionless 2D Toda hierarchy of partial 
differential equations, so that coefficients of the expansion satisfy an infinite set of non- 
linear relations. Also, we have seen that the solution of the Toda hierarchy is selected by a 
non-trivial string equation v _1 } = 1. The following explicit results have been found: 

o First, one can compute some infinite series of coefficients of the expansion using 
these differential equations. The example worked out in section 4 is the set of coefficients 
of the form 9\9\\ if x = 0i0i, then £ = xdF/dx was found to satisfy the equation 16£ 3 + 
8£ 2 + (1 - 36x)£ + x(27x - 1) = 0, that is 

F = Mi + i(0i0i) 2 + ^(0i0i) 3 + 6(0!0!) 4 + ^(0i0i) 5 + • • • (6.1) 

More generally, one can compute generating series of coefficients involving only a 
finite number of distinct 9 q and 9 q as solutions of algebraic equations; a more complicated 
example was given without proof, cf Eq. ( |4. 13| ) . In principle, this allows to compute 
any given coefficient required, although it does not give a closed expression for the whole 
function of the infinite set of 9 q and 9 q . 

o Secondly, one can obtain all the coefficients (as functions of the 9 q ) of the first 
few orders of the expansion in products of 9 q : we rewrite their generating functions here 
(Eqs. O, Q and fl03|) ): 

^-V(a)F = b(a)-~ (6.2.1) 
da a 

V( 0l )VMF = logM^M (6.2.2) 

1/ai - l/a 2 

, . , . , . d&(oi)/doi db(a 2 )/da 2 

V(ai)V(a2)V(a 3 )F = — — r — — - — r — -- + 



(6(ai)-6(a 2 ))(6(oi)-6(a3)) ip(a 2 ) - b{a x )){b{a 2 ) - 6(a 3 )) 

, db(a 3 )/da 3 . , 

+ (6(o3)-6(oi))(6(o3)-6(aa)) 1 ' " ' 

where in the rescaled variables, V(a) = ^2 q >iCL q -^~- The expansion is given in terms 
of the function 6(a) = 1/a + ^ g>1 m q a q ~ l . The free cumulants m q can be reexpanded 
themselves in powers of the 9 q . fh q = 9 q — |r ij>i 9i9j + ^2 i,j,k>i 9~i9j9~k + •■■■> 

i + j = q ' i+j + k = q 

so that we find for example 
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- 2 ^Al^ 

q,ri,r 2 >l 
ri+r 2 =g 

- 2 QqiQq^r 

r,qi,q2>l 
qi+qz=r 

+ ^ (?i + 92-nun(gi,52,ri,r 2 ) + l)^ 1 ^ 3 6'ri6'r2 

gi,g 2 ,T-i,r 2 >l 
<2i+<?2=T-i+r 2 

+ g ^ (? + ^)0q9 ri 9r 2 9 r3 

q,ri,r 2 ,r 3 >l 
g=n+r 2 +r 3 

+ i (r + l)9 qi 9 q2 9j r 

r,qi,qi,qs>l 
r=qi+q2+qz 

+ ■■■ (6.3) 
with the appropriate symmetry of interchange of the 9 q and 9 q . 
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